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Abstract 

We verify the functional form of the asymptotics of the spin-spin equal -time 
correlation function for the XX-chain, predicted by the hypothesis of conformal 
invariance at large distances and by the bosonization procedure. We point out 
that bosonization pocedure also predicts the functional form of the correlators 
for the chains of finite length. We found the exact expression for the spin-spin 
equal-time correlator on finite lattice. We find the excellent agreement of the 
exact correlator with the prediction given by the leading asymptotics result up to 
the very small distances. We also establish the correspondence between the value 
of the constant before the asymptotics for the XX-chain with the expression for 
this constant proposed by Lukyanov and Zamolodchikov. We also evaluate the 
constant corresponding to the subleading term in the asymptotics in a way which 
is different from the previous studies. 



It is well known that the hypothesis of conformal invariance at large distances predicts 
both the critical exponents and the functional form of the correlation functions for 
massless one - dimensional systems on the circle of finite length L (for example see 
In this approach the conformal mapping w = ^\n{z) of the infinite plane z on the 
strip of finite width w = t + ix, x ^ {^,L), is employed. It is known that in the case 
periodic boundary conditions (for the initial spin operators) the low-energy theory - is 
the Conformal Field Theory with cental charge c=l (which can be computed for XXZ 
chain for an arbitrary boundary conditions), which is equivalent to the so called Gaussian 
model of CFT with c=l, with known spectrum of primary operators. The only parameter 
is the compactification radius (see below). I would like to stress also that application of 
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CFT for calculation of e.g. critical indices is the hypothesis in a sense that it is not proved 
in a mathematically rigorous way. Thus any independent verification of the conformal 
field theory predictions in the exactly solvable models is of interest. Alternatively, one 
can use the bosonization procedure for the low-enegry effective theory for the systems 
in the gapless regime (Luttinger liquid) in its different versions 0,0 to predict the 
critical exponents for various one- dimensional systems (for the first calculation of the 
critical index for the XXZ -spin chain see ref. ^). In our opinion it is obscure in the 
literature that the bosonization also predicts the functional form of the correlators in 
the form ~ cos(2pi;'x)((l/L)sin(7rx/L))", where L -is length of the system, a is the 
critical index and - is the Fermi-momentum, at the sufficiently large distances. So it 
is worth pointing out how this functional form appears in the bosonization framework. 
This is done in the first part of the paper. Next, it is interesting to verify this functional 
form in the model where the correlations functions can be calculated exactly. In the 
present paper we study the spin-spin equal-time correlator in the XX - spin chain (the 
hamiltonian is if = \ Yl!i=\{p^ + c^fcrf+i)) the anisotropy parameter A = 0). It is 
worth mentioning that although the XX-chain can be solved with the help of mapping to 
the free-fermion system via the Jordan- Wigner transformation the system is essentially 
the the hard -core bosons on the lattice (at the half-filling) and is not equivalent to the 
free fermions which manifests for example in the completely different correlators in the 
two models. We show in the present letter that the spin-spin correlator in the XX- model 
can be calculated exactly on a chain of finite length at any distance x of order of the chain 
length (see eq.(^ below). We show that the functional form of the correlator given by 
d^) coincides with that given by the conformal fild theory (bosonization) . Note also that 
the knowledge of the functional form of the correlators can be usefuU in the numerical 
study of correlators and extracting the critical exponents in different models. In Section 
1 we show how the usual bosonization procedure leads to the functional form predicted 
by the conformal invariance. In Section 2 we brifly review the well known calculation 
of the spin -spin correlator in the XX -model [^], and present the expression for 
the correlator at the distances x ~ L. Finally in Section 3 we compare the numerical 
constant in the asymptotic of the correlator given in [Q with the expression for this 
constant proposed recently in ref. and evaluate the coefficient corresponding to the 
subleading term in the correlator 0, in a new way. 

1. Bosonization. 

Consider the effective low -energy Hamiltonian which build up from the fermionic 
operators (a/c, Cfc, k = 2TTn/L, n ^ Z, L- is the length of the chain) corresponding to 
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the excitations around the right and the left Fermi- points and consists of the kinetic 
energy term and the interaction term H = T + V with the couphng constant A: 

H = Y1 Kot(^k - ctck) + A/L ^ alak+qct,Ck'-q. (1) 

k k,k',q 

Defining the operators 

k k 

where \k\, \k+p\ < A, where A is some cut-off energy, which for the states with the filled 
Dirac sea have the following commutational relations 

[pii-p); piip)] = ^^p,p' [p2ip); p2i-p)] = ^K,p'^ 

one can represent the Hamiltoniamn in the following form: 



2n 2tt 

H = -j-Yl {P1{P)P1{-P) + P2{-P)P2{P)) + ^Y.-r (P1(P)P2(-P) + P1(-P)P2(P)) • 
^ p>0 p>0 ^ 

First, for the case of free Dirac fermions A = the system is equivalent to the free scalar 
Bose field with the Lagrangian L = ^{{dt4>y — {dx(pY) 



cPix,t) = ^Y.^{bte"''' + bke- 
VL 2J\k\ ^ 



ikx 



6kk' and kx 



where bk,b^ are Bose creation and annihilation operators hki]h^ 
\k\t — kx. The free Hamiltonian takes the following form: H = J2p\p\bpbp, where 
p = 27m/L and the equal time commutation relations are [dt4>{x); (j){y)] = —i6{x — y). 
The identification of these operators with the operators pi.2{p) defined above is 



> 0) = p,{p)^27r/L\p\ b;{p < 0) = p2{p)^2n/L\p\ 

for the free fermion case. 

To evaluate the correlators in the system of finite length and make the connection 
with the conformal field theory predictions, one can proceed as follows. First one can 
define the lattice fields ni 2(0;) with the help of the Fourier transform as 

PiM =Y.^'^''^h2ix)^ "1,2(3;) = ]rJ2e-''"'pi^2{p) 

X ^ p 

This fields have a physical meaning of the local number of the fermions above the Fermi 
level at the right and the left Fermi points. In terms of this fields the Hamiltonian has 
the following form: 

H = 2t:Y^ +^2(3;)) + Xni{x)n2{x) 
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Considering the average distribution of the number of extra particles we obtain the 
ground state energy in the form 

= ^ Q ((AiVi)2 + {AN2Y) + XAN,AN2^ 

where AA^i 2 - are the additional numbers of particles at the two Fermi points. One can 
also rewrite the ground state energy in the sector with the total number of particles and 
the momentum AA^ = AA^i + AA^2, AQ = ANi - AA^2 in such a way that the total 
Hamiltonian takes the following form (this for was first proposed in ref.[^]): 

H = u{X) E IpHK + ^u(A) [am' + (1/0 (AQ)^] , (2) 

where the parameters u{X) = (1 — A^)^/^ and ^ = ((1 + A)/(l — A))^/^. Calculation of 
finite - size corrections to the energy of the ground state for the XXZ- spin chain (see 
for example [|l^) leads to the expression (|) and allows one to obtain the parameter ^ 
which leads to the predictions of critical indices according to the conformal field theory. 
The calculation gives the value ^ = 2{7i — r])/^, where the parameter 77 is connected with 
the anisotropy parameter of the XXZ - chain as A = cos{ri). Next one establishes the 
commutational relations for the fields ni^2{x)- We have: 

To take the continuum limit it is sufficient to extend the sum over the momentum to 
p G (±C)o) (note that initially it was implied the sum in the limits p G {±tt)). Then we 
obtain 

[ni{xy,n,{y)] ^ J dppe'^^^-y^ = -1-6' {x-y) 

Then intoducing the new variables 

^1,2(2;) = y27r ni,2(x), 
we have the following commutational relations and the density of the Hamiltonian 

H = \: (ni(x)ni(x) + h2{x)h2{x)) + Ani(x)n2(x). (3) 

We also have the following commutational relations \ji\[x)] n\[y)\ = S'{x — y). We then 
have the following conjugated field and the momenta: 

7r(x) = -^{ni{x) - n2{x))] d^(t){x) = -^{ni{x) + n2{x)) 
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= N{x) = Ni{x) + N2ix), N,^2ix) = / dy 

Jo 

In terms of these variables the Hamiltonian takes the following form: 

H = ^£dx ((1 - \)n^ix) + (1 + \){d(P{x)f) . 
The Hamiltonian density can also be represented in the following form: 
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H = -u{X) + ad<l>{x)y\ = -u{X) [n\x) + {d<P{x)f\ , (4) 

where 

7r(x) = V^7r(a;), <P{x) = {1/ foH^)- (5) 
The last equation (|^) is nothing else as the canonical transformation, which is equivalent 
to the Bogoliubov transformation for the original operators pi.2{p)- Next to establish 
the expressions for Fermions one should use the commutational relations [a"^(x); pi{p)] = 
—e^P^a'^{x) and the same for c'^{x). Note that these last relations were obtained using 
the expression with original lattice fermions: pi{p) = Yl,y ^-^^^ {y)a{y) . In this way we 
obtain the following expressions for fermionic operators: 

a ' (X) = i^i exp I ^ ^ PiMg-*px j _ {-i2TtNi{x)) , (6) 




c+(x) = exp - — 5J ^IfMZg-^P- = K2 exp {i2'kN2{x)) , 

where the fields Ni 2{x) differ by the normalization \/2tx from the fields iVi 2(0;) and K12 
are the Klein factors - the operators which creates the single particle at the right (left) 
Fermi -points (we omit here the usual exponential suppression e~°'^'/^ in the exponent 
and and the constant factor in front of the exponent l/VSvra which in the limit a ^ 
leads to the correct anticommutational relations [^]). Note that the above expressions 
(13) are equivalent to the known "field-theoretical" bosonization formulas 0: 



a (c^)(x) = Ai,2 e ^ ' 



Now let us apply the above formulas to the specific case of the XXZ - spin chain. 
Using the Jordan- Wigner transformation = a^exp{i7rN{x)), where stands for the 
"original" lattice fermionic operator, and performing the obvious substitutions N{x) 
x/2 + Ni{x) + N2{x) and a+ e*^^^a+(x) + e~*''^^c+(x), pp = vr/2, we obtain after 
the canonical transformation (^ the expression for the spin operator which determines 
the leading term in the asympotics of correlator for the XXZ -chain: 

at ~ (-l)^exp i^-^nfaNl{x) - M^))) , (7) 
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where iVi_2(a;) - are corresponds to the free fields vr(x), 0(x), obtained after the trans- 
formation To these operators correspond the new operators pifi{p) and the new 
fermionic operators (quasiparticles). Analogously the term responsible for the sublead- 
ing asymptotics has the form 




Averaging the product of exponents in bosonic operators for the expression (|^) and using 
the properties of pi,2(p), (pi(-p)pi(p)) = ^d{p) and (p2(p)p2(-p)) = 2^6'(p), we get for 
the correlation function G{x) = {0\a^_^^o-^ \0) the following sum in the exponent: 

C expl - y le^"(2-x/L) ^ \ 

where C - is some constant. Then using the formula ^-2" = — ln(l — z) and 

substituting the value ^ = 2(7r — ?7)/7r ^ 1 we obtain the following expression for the 
XX - chain: 

G{x) = Co , a = ^^ = 1/2 ix»l). (8) 

(Lsin(f)) TT 

Thus, although bosonization, which deals with the low-energy effective theory, is not able 
to predict the constant before the asymptotics, the critical exponent and the functional 
form are predicted in accordance with conformal field theory. 

2. Asymptotics of correlators in XX -spin chain. 

Let us briefiy review the exact calculation of the spin-spin equal-time correlation 
function (density matrix) for the XX- spin chain on finite lattice of the length L: 

^(x) = (0|(T+,ar|0). 

Using the Jordan- Wigner transformation relating spin operators to the Fermi operators 
{af, ai) (7+ = exp{iTTj2i<x^i)(^ty correlation function G{x) can be represented as 
the following average for over the free-Fermion ground state: 

G{x) = (0|a+e*-^(^)ao|0), 

where N{x) = J2i=i n-i. Introducing the operators, anticommuting at different sites, 

Ai = at + ai, Bi = at-ai, AiBi = e'^''^ 
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where Ui = afai- is the fermion occupation number, with the following correlators with 
respect to the free-fermion vacuum, 

{0\BiAj\0) = 2Goit-j), {0\AiAj\0) = 0, (0|5,5,|0) = 0, 

where the free-fermion Green function on finite chain Gq{x) is 

n f \ /nl + ln\ sin(7rx/2) 
Go(x) = {0\at.am = J—^y 

one obtains the following expression for the bosonic correlator: 

G{x) = ^{0\Bo{A,B,){A2B2) . . . (A,_i5,._i)A,.|0). 

Note that we assume the periodic boundary conditions for the initial spin operators, so 
that strictly speaking, the above formulas are valid only for the case when L - is even (the 
ground state is not degenerate) and M = L/2 - is an odd integer (so that L/4 - is not 
an integer). In fact it is easily seen that for the boundary term ~ exp(z7r(M — 1)) = +1 
when M - is odd, the momentums of fermions - are integer (not half-integer), their 
configuration is symmetric around zero and the free - fermionic Green function Go{x) - 
is given exactly by the above formula. It was also mentioned in ref. |0 that the periodic 
boundary conditions for spin operators ("a-cyclic" in this paper) leads to the same 
results as for the periodic boundary conditions for the fermionic operators ("c-cyclic") 
for the case when L/4 - is not an integer. Using Wick's theorem we obtain the following 
determinant oi x x x matrix: 

G{x) = detij{2Go{i - j - 1)), i,j = l,...x. 

Due to the form of this matrix {Go{l) = for even /) this determinant can be simplified 
and the following formulas are obtained: 

G{x) = ^{RnY, {x = 2N), G{x) = {x = 2N + 1), 

where we denote hj Rn the following determinant of the x A^- matrix: 

Rn = det,,((-l)^-^Go(2« - 2j - 1)), j = 1, . . . iV, 

where Go{x) is the same Green function of free fermions as above. Since Rn - is the 
Cauchy determinant one can obtain the following expression for it on the finite chain: 

\7iJ ^},\sm{7c{2k + l)/L)sm{n{2k-l)/L)J ' 
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This formula is the main result of the present paper. Note that we obtained the ex- 
act expression for the correlator on the finite lattice. As it was mentioned above the 
expression (j^) for the correlator was obtained for M - odd. As for the case M - even, 
it is clear on general grounds that in this case the correlator G{x) will be modified by 
the terms of order 1/L at any distance x (that can be proved rigorously, however the 
proof will not be considered in the present letter). From the expression (^) it is easy 
to obtain the correlator in the thermodynamic limit (L oo) which is given by the 
similar product. We consider the analitical evaluation of this product for large in the 
next Section. For the finite chain it is easy to evaluate (^) numerically and compare 
the result with the asymptotic (Q) at x >> 1 and x ~ L. One finds numerically that 
already at the sufficiently small distances x, where the contribution of the subleading 
term ~ becomes negligible, the correction to the asymptotic formula behaves 

like ~ 1/L (at the very large L ~ 10^ it is difficult to evaluate the product at x ~ L due 
to the numerical reasons). We find numerically that the exact correlator coincides with 
the correlator given by (|^) with very high accuracy up to the very small distances x ~ 1. 
One can also consider the limit L — oo in (|^) and using the exact value of G'(x) at x ~ 1 
estimate the value of the subleading term which is predicted to be of order ~ l/x"^^/"" 
for the XXZ - chain. The coefficient for this term appears to be very small: at x ~ 1 
(see Section 3, equation (p!3D) so that the correlator for the infinite chain can be seen in 
a very good agreement with the asymptotic formula (^ even at the very small distances 
X ~ 1. The contribution of the next terms is difficult to obtain numerically (presumably 
due to their smallness). The physical reason for the smallness of the subleading term is 
not clear at the moment. 



3. Numerical constant. 



Here we calculate the constant Cq in front of the asymptotics of the density matrix in 
the XX -spin chain and derive the representation proposed in ref. for this case. Using 
the expressions in the thermodinamic limit we calculate the product: 

^"-l^j ll\(2k+l)(2k-l)) ■ 

Considering the logarithm of i? at we obtain the following sum: 

ln(i?^)=Arin(2/vr) + X:-(2) Jli^ - k)^^ 
Next, we use the following formula for the finite sum YJk=i l/^"* = C(s) — Ci^^n + 1) 



8 



(where ({s,n) -is generalized Riemanian zeta- function) |T2[: 

where ijj^'^^z) is the ra-th derivative of ^/'-function ^p{z) = V' (z)/T(z), which leads to the 
following expression for ln(i?7v): 

-'VEi(i) ^*-'->W-gi(i) (10) 

since the sum 

cancells the term coming from the factor {2/tt)^ . The last formula can be obtained using 
= (— l)"+^n!C(n + l) and rewriting the last sum using the definition of function 
as — I^fcLi ^ 1/4/c^) = — A^ln(2/7r), which can be seen from the infinite product 
representation sin(x)/x = nfcli(l — tt'^x'^ /k"^). Thus we are left with three terms in 
eg. (p!OD . To find their asymptotic behaviour at large N up to the terms of order ~ l/N"^, 
we use the asymptotics 



1 ^ B. 



i^{z) = Hz) - - - E 



2n 



2z 2nz2" ■ 



\z\ CX), 



where - are Bernulli numbers {B2 = 1/6). One can see that the sum of the second 
and the third terms in eq.(|IO|) is 

while the first term equals 

where 7 = = 0.5772 ... - is the Euler's constant. To estimate this sum one can 

use the integral representation: 



^{z) = dt{ 



1-e 
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One also readily obtains the integral representation for il)^"^^ '^\^) which leads to the 
following expression for the first term: 




Thus combining all terms after the simple algebra, using the integration by parts and 
taking into account the constant term —1/4 in eq.(|l^), we finally obtain the result: 

. _i,„(iV) + i f f (e-« - ^) - (12) 

where the omitted terms are of order ~ 1/N'^. One can see that the above expression 
coincides with the formula proposed in ref. in the case of the XX - chain. The last 
term in ( [T^) gives the following coefficient for the next-to-leading asymptotics for the 
correlator: 



where the constant Co is defined in (||). The value of the constant corresponing to the 
subleading term in eq. (|T3[) was first obtained by McCoy [^] using the method with the 
help of the asympotics of the Barnes G- function |Tl|], defined by G{z + 1) = V{z)G{z), 
G{1) = 1. The product Rn can be represented as 



G(A^+ l/2)G(iV + 3/2)' 
Using the asymtotics of the function G{N) at large A^, 

G{N) = — -\nA- -(ln27r)A^ + (-N^ - —] InN - -N^ + of— 



12 2' ' V2 12/ 4 VA^2 

where A is the Glaishers constant (see below), one can obtain the result ([T3| ) with the 
constant (Cq)^^^ = n^^^2^^^'^e^^^A^'^ . Note that using the asymptotics of ip{z) one can 
find analytically the higher terms in the asymptotics of the correlator. One could also 
represent the first term in eq.(|lOD as 

which also allows a numerical estimation of the constant. One can evaluate the integral 
in (|12D to get the asymptotic 



10 



which is equivalent to the estimate 

= (iV-V4) (2Vi2ei/4^-3j ^ 

where A- is the Glaisher constant 

A = ei/i2-C'(-i) = 1.282427.... 
This result agrees with the result obtained by Wu using the expression of the product 



through the Barnes G- functions |TT[. Finally, an alternative way to obtain the integral in 



eq.(|T2l) is to use the relations for T - function T{n+1) = n\, T{n+l/2) = 7T^/^2~'^{2n-l)\\ 
to represent in the form: 



^ tiT{k + l/2nk-l/2y 



Then one can use the known integral representation for the sum of the type Ef=ilnr(fc + 
b) (see.|jl^), and extract the large asymptotics of the resulting expression. Thus for 
the constant before the asymptotic (H) the value Co/2-\/7r = 0.147088 ... is obtained. 

Conclusion. 

In conclusion, we verified the functional form of the asymptotics of the spin-spin 
equal -time correlation function for the XX-chain predicted by different methods. We 
find the excellent agreement of the exact correlator with the prediction given by the 
leading asymptotics result up to the very small distances. We have also estimated the 
coefficient corresponding to the subleading correction p , in a different way and found 
the expression for the leading term in agreement with ref.0. It is worth mentioning that 
the correlators for the XY- spin chain and the Ising- like spin chain have been studied 
at the non-zero temperature and in the time - dependent case in different limits in a 
number of papers (for example, see ref. |]13|). Let us mention that the exact correlator in 
the XX- model (at least in the infinite L limit) can be obtained in the framework of the 
general approach based on the Algebraic Bethe Ansatz method (see for example [|T^) 
as well as in the approach based on the formfactors of the model [jl5[ which is beyond 



the scope of the present letter. Let us mention also that studying of the formfactors 
with the small energy of the excited states can be the basis of the determination of 
the constant Cq. We postpone this question for a separate publication. It is possible 
that the formfactor approach can be usefuU for the calculation of the constant Cq in the 
general case of the XXZ - spin model. 
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